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Abstract 

Experiments observing the liquid surface in a vertically oscillating container have indicated 
that modeling the dynamics of such systems require maps that admit states at infinity. In 
this paper we investigate the bifurcations in such a map. We show that though such maps 
in general fall in the category of piecewise smooth maps, the mechanisms of bifurcations 
are quite different from those in other piecewise smooth maps. We obtain the conditions 
of occurrence of infinite states, and show that periodic orbits containing such states are 
superstable. We observe period-adding cascade in this system, and obtain the scaling law 
of the successive periodic windows. 
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1 Introduction 



Recently a lot of research attention has been directed toward the dynamics of piece- 
wise smooth maps (PWS), because they represent a large number of systems of 
practical interest including switching electrical circuits and impacting mechanical 
systems. In such systems the discrete phase space is divided into compartments 
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within which the map is smooth, and the compartments are separated by border- 
lines at which the map is not differentiable. A one-dimensional piecewise smooth 
map has the general form 

g(x n ,fi), for x n <\ 
Xn+i = t{x n ) = \ (1) 
h(x n ,/j),for x n >\ 



where // is the bifurcation parameter and the compartments are separated by the 
borderline value A. In such a map, there is the possibility that a fixed point may 
collide with the border with the change of a system parameter. When that happens, 
there is a sudden change in the stability of the fixed point, leading to a new type of 
nonlinear phenomenon called border collision bifurcation [1,2]. It has been shown 
that such border collisions may lead to atypical bifurcation phenomena like transi- 
tion from period-2 to period-3 or a sudden onset of chaos without undergoing the 
usual period doubling cascade [3]. 

A few different forms of such maps have been investigated to date: 

(1) The map / is continuous, not differentiable at A, and both dg/dx and dh/dx 
are finite. Such maps represent a class of switching circuits, and have been 
investigated in detail [2,4]. 

(2) The map / is continuous, not differentiable at A, and there is a square root sin- 
gularity (i.e., dh/dx is infinite) at one side of the border. Such maps represent 
the impact oscillator [5,6,1]. 

(3) The map / is discontinuous at A, the derivative df/dx is also discontinuous 
at A, but the value of the derivative at both sides of the border are finite. Such 
maps represent a class of electronic circuits [7] including the Colpitts oscilla- 
tor [8], and the bifurcation theory for such maps has been developed recently 
[9]. 

In 1997 an experiment was reported, where the oscillations in the surface of a liquid 
held in a vertically oscillating container were observed using a laser probe, and it 
was found that under some conditions narrow jets are ejected from the center of 
the surface. It was shown that representation of this system required a map with 
not only slope singularity but also magnitude singularity at the border [10]. The 
proposed map had the form 

for x n < A (2) 
for x n > A (3) 

where a, (3, p and A are constants and 7 is the bifurcation parameter. The graph 
of the map is schematically shown in Fig. 1. In this paper we investigate the bifur- 
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cation phenomena in a map of the above form - for which no theory is currently 
available. 




X x n ► 

Fig. 1. The graph of the map given by (2) and (3). 



In this map, the vertical line x = A forms an asymptote for the equations (2) and (3), 
and a singularity occurs at this value. This asymptotic behavior occurs due to the 
geometric considerations in the waves — that the waveheight/wavelength cannot 
exceed some ratio before the wave becomes self-intersecting. While obtaining the 
bifurcation diagram, if x takes a value close to A, the value of x at the next iterate is 
very high. The program has to account for this possibility. The bifurcation diagram 
obtained this way is presented in Fig. 2. As both domain and range of the piecewise 
smooth map are represented by the half-open set [0, oo), Fig. 2 and all subsequent 
bifurcation diagrams have been truncated. 
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Fig. 2. Bifurcation diagram with a = 0.04, (3 = 2.8, p = 0.14, A = 5.76, and 7 as the 
variable parameter. The diagram is truncated above x n = 10. 
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A few features are noticeable in the bifurcation diagram. First, there is a succession 
of periodic windows with the following properties: 

(1) Within each periodic window, as the bifurcation parameter 7 is reduced, a 
period doubling cascade terminates in chaos. 

(2) If n is the periodicity of the base orbit (lowest period orbit) at the highest 
parameter value in a perodic window and m is the periodicity of the base orbit 
at the highest parameter value in the next window, then m = n + 1, i.e., there 
is a period adding cascade as the parameter is reduced. 

(3) The width of the periodic window reduces monotonically as the period adding 
cascade progresses. 

In earlier studies, period adding cascades were observed in the study of piecewise 
smooth maps of finite magnitude and finite slope [2,4] and those with square root 
singularity of derivatives [5,6,1]. In maps of the former type it has been found that 
each periodic window originates at a border collsion. But for the system under 
the present study, the bifurcation diagram plotted for a larger range of x (Fig. 3) 
demonstrates that each periodic window does not emerge due to border collision. 
These occur at saddle-node bifurcations. 

In systems with square root singularity it was observed that a period-n orbit always 
has n — 1 points on one side of the border and one point on the other side. Such 
orbits were called maximal. Though from Fig. 2 this may seem true also for the 
system presently under consideration, a closer scrutiny of Fig. 3 shows that the or- 
bits are not maximal at the points of emergence of periodic windows. They become 
maximal as the periodic point crosses the borderline value of A. 
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Fig. 3. Bifurcation diagram with horizontal lines drawn at y = j3 and y = A. 
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2 Border-collision bifurcations 



In view of the differences, therefore, the theory developed for these other types of 
piecewise smooth maps cannot be applied for the map with magnitude singularity 
presently under consideration. 

For this system, a few questions demand fresh answers. For which parameter values 
does the state reach infinite value? It is easy to see that in the chaotic windows, 
because of the ergodic nature of the orbit, iterates come arbitrarily close to the 
value of A, and then the next iterate must assume infinite value. Therefore, even 
though a bifurcation diagram cannot be drawn for such high range of the variable, 
there must be points at infinity for all parameter values in all the chaotic windows. 
This is supported by the observation in [10] where such states were called "ejecting 
states." 

This system reveals another interesting aspect: infinity is reached even within peri- 
odic windows — at the points of border collision. 

The points of border-collision in the periodic windows can be determined theoreti- 
cally. At the point of border collision, one point of the periodic orbit must have the 
value A. Since A i— > oo i— > (3, for a period three orbit the three points are (3, A, oo. 
To find the value of the parameter 7 at the crossover point we set 

^ + (^ = * (4) 



This yields a value 7 = 2.052577 for which the period-3 orbit contains a point at 
infinity. Fig. 4 gives a closeup of the period-3 window showing this point. 

At the point where the period-4 orbit touches infinity, three of the four points will 
be A, 00, and (3. The fourth point obtained by substituting x n = (3 in (2) is 

113 + 



Since this point must map to A, we obtain the parameter value 7 from (2) as 



(/?-A) 2i (7/3+^-A)^ 



7{7^ + + :r )2 L = a (5) 



This yields the following 4-th degree equation in 7 

a 7 4 + ai7 3 + a 2 / y 2 + a 3 7 + a 4 = (6) 
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where 



a = (3 3 



Fig. 4. Close up of the period-3 window. 
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Solving (6) using (3 = 2.8, A = 5.76, a = 0.040 yields 7 = 1.425510183648710. 
Out of the other three solutions, one is negative and the other two complex conju- 
gate — which are not possible values of the parameter. 



In a similar manner, all the periodic orbits with a point at infinity and their corre- 
sponding parameter values can be obtained. The only requirement for such border- 
collision points to exist is (3 < A. 
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3 Stability of periodic orbits 



The rate of expansion/contraction at any periodic point is dependent on the deriva- 
tive of the map function 



= 1+1 - 7 for x n < A, (7) 



dx n (x n - A) 2 (x n - A) 3 

dxn+i _ p 2px n 

dx n [x n — A) 2 {x n — A) 3 



for x n > A. (8) 



Since border collision points are related to the singularity condition, and since from 
(7) and (8) we find that 

dx n +i dx n+ \ 
hm — > oo and hm — > — oo, 



x n >x dXyi x n ^a^~ dx 



n 



one may tend to believe that periodic orbits containing a point at infinity cannot be 
stable. However, since 

lim rf*n+l ^ 



this need not be true, and so it is interesting to work out the stability of the orbits 
containing infinity. 

Such an orbit must have two points with values A and oo. The other points in the 
orbit have finite slope. Therefore the derivative of nth iterate map for that condition 
must be a finite number times the limiting values of 

,. , ot ax n \ ( p 2px n+ i \ 
lim 7 + t TTf - 7 rig x lim 7 rrf ~ 1 TTa ( 9 ) 

Xn ^\- {x n -xy {x n -\yj \x n+1 ^co {x n+ i - \y (x n+ i - \y j 



and 



Hm P 2px ra \ ( p 2px n+1 

; n ^A+(x n -A) 2 (X n -A) 3 / \x n+1 ~^oo (x n+ i - A) 2 (Xn+i-A) 3 , 



For obtaining the limiting values of (9) and (10), we transform these into a limit 
of a single expression. Using the symbolic computation facility of Matlab, we 
substitute the expression for x n+1 , i.e., (2) or (3) depending on the position of x n . 
For simplicity and compactness these expressions are denoted by F~ and F + as 
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x n — > A or x n — > A + . The limit expression becomes 

Mm F- (11) 

lim F + (12) 

We find that in both cases limit exists and is zero. Therefore we conclude that the 
orbit including the point at infinity is superstable. It should be noted that the value 
of limit is independent of (3, A, a, p and hence the result is a general one. 

The result was numerically confirmed by choosing a border collision point and then 
by calculating the product of the derivatives. It was observed that the product could 
be made arbitrarily small by choosing appropriate value of the variable x. 



4 Bifurcation behavior for (3 > A 
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Fig. 5. Bifurcation diagram with respect to (3 for (3 > A. 

In the earlier sections the bifurcation behavior for the case (3 < A has been in- 
vestigated. We now consider the situation where (3 > A, and study the bifurcation 
phenomena as (3 is varied. The value of 7 only changes the graph quantitatively 
and hence it is not considered to be the parameter. Moreover, p is assumed to be 
positive. 



The behavior can be explained in the light of superstability of periodic orbits when 
border-collision occurs. When (3 = A the points A, 00 form a superstable 2-period 
orbit. As (3 is increased, the product of the derivatives approaches unity. At a point 
the period-2 orbit becomes unstable giving way to a stable period- 1 orbit. This is 
a standard period-doubling bifurcation when (3 is reduced. The nature of the map 
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(Fig. 1) also suggests that period- 1 orbit should be stable at very high values of (3. 
This behavior is shown in the bifurcation diagram of Fig. 5. 

However, if the initial iterate starts on the left-hand side of the map and remains 
forever on the left-hand side, then the behavior as shown in Fig. 5 is not exhibited. If 
< 7 + < 1 then the origin forms a stable fixed point and iterates starting on 
the left hand side of the map converge to the origin. 



5 Scaling in the period-adding cascade 

Figure 2 reveals that successive periodic windows have monotonically decreasing 
width, which suggests the existence of a Feigenbaum-type ratio. To check this, we 
obtained the parameter values at the points of saddle-node bifurcation where the 
periodic orbits come into existence. The ratio of the widths of successive periodic 
windows were taken. Fig. 6 shows the graph of the ratios of the widths versus the 
index. 




1 1 1 1 1 1 1 1 1 

0. 2 4 6 8 10 12 14 

Index 



Fig. 6. The ratios of the widths between the successive appearance of periodic windows. * 
indicates the numerically determined values of the parameter ratios, and o indicates those 
obtained from (13). 

Using standard curve-fitting technique, we fit the data into the curve 

y = ax b + c. (13) 

The obtained value of c was 1.091 ± .033. Hence the ratio of successive widths will 
tend to c. 
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6 Conclusions 



In this work a piecewise smooth map was considered that has magnitude as well 
as slope singularity at a borderline value. It was found that the singularity condi- 
tion, i.e., the state assuming infinte value, occurs both in chaotic mode as well as 
in periodic mode. We found that periodic orbits containing a point at infinity oc- 
cur are the points of border collision, at which these orbits are superstable. The 
system exhibits period-adding cascade with deminishing width of successive peri- 
odic windows, and the ratio converges to a number 1.091 ± .033. There are several 
border-collisions within each periodic window. 
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